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Abstract 

In this paper a procedure for large-eddy simulation (LES) has been devised for fluid and magne- 



tohydrodynamic turbulence in Fourier space using the renormalized parameters. The parameters 
calculated using field theory have been taken from recent papers by Verma 



,|2j. We have carried 

out LES on 64 3 grid. These results match quite well with direct numerical simulations of 128 3 . We 
show that proper choice of parameter is necessary in LES. 

PACS numbers: 47.27.Eq, 47.65. +a, ll.10.Gh 
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Turbulence is one of the most difficult and unsolved problems of classical physics. To 
probe the complex dynamics of turbulence, one often resorts to computer experiments, 
known as Direct Numerical Simulations (DNS). Since multiple scales are involved in tur- 
bulence, DNS of turbulence is a very expensive task in terms of both computer time and 
memory, even in modern computers. For example, a pseudo-spectral simulation by Gotoh [^] 
on 1024 3 grid using vector parallel Fujitsu VPP5000/56 with 32 processors took 500 hours of 
computer time, and required 8 Gigabytes of memory per processor. To reduce the required 
computer time and memory space, an ingenious technique called large-eddy simulation (LES) 
has been developed (see review article by Metais |4( and references therein). 

Basic idea of LES is to resolve only the large scales of turbulent flow. The effect of 
smaller scale interactions are modeled appropriately using the existing theories. In turbu- 
lence, Fourier modes of different scales interact with each other. Kolmogorov provided an 
important model of turbulence in which the interactions effectively yield a constant energy 
flux from large scales to intermediate scales, and then to small scales. When we observe 
Fourier modes up to certain length scale I in the intermediate range, the modes with scales 
less than I act as a sink of energy. According to Kolmogorov' s theory, the amount of sink 
should be equal to the energy flux. In LES, the large scales up to I are resolved by using 
eddy viscosity at cutoff scale /, where energy is drained. Analysis of turbulence using renor- 
malization groups (RG) shows that the above modeling is possible. LES uses this idea to 
analyze large-scale dynamics of turbulence. 

Renormalization Group (RG) is a popular tool used by physicists to solve problems with 
multiple scales. Since turbulence involves multiple scales, RG has been applied successfully 
to turbulence . In Wilson's Fourier space RG scheme, Fourier space is divided into 

many shells. The nonlinear interactions among various shells are computed using first-order 
perturbation theory, that yields an effective viscosity, called renormalized or eddy viscosity, 
at any scale. The renormalized viscosity is found to be wavenumber (k) dependent. McComb 
and Watt computed the renormalized viscosity using 'self-consistent' RG procedure. 
When the cutoff wavenumber kc is in the inertial range, the renormalized viscosity is given 
by 

Mkc) = (iO^n^VV (l) 

where II is the energy flux, K is Kolmogorov's constant, and v* a parameter. McComb and 

n 

Watt 8j found ~ 0.50 and K ~ 1.62. Verma []J also computed the above quantities using 



a refined technique and found v* ~ .38 and K fa 1.6. 



Zhou and Vahala |9|, [l0| developed an alternative recursive-renormalization-group theory 
for turbulence modeling. In their calculation they find backscatter of energy from small 
scales to large scales, and a cusp in renormalized viscosity near kc- These features are 
attributed to triple correlations, which has not been accounted for in McComb and Watt's 
calculations. Recently Schilling and Zhou 11] have addressed the above problem using eddy- 
damped quasinormal Markovian (EDQNM) closure model. In the current paper we neglect 
backscatter. 

McComb [7] had proposed that the renormalized viscosity v(kc) could be used as effective 
viscosity for LES, however, this calculation had not been done till date. Earlier, the spectral 
eddy viscosity u t (k\kc) has been used for LES in EDQNM formalism (see [l2j). In this 
scheme, 



Kc) = 0A4lK- 3/2 



E{k c ^' 2 



f(k/k c ) (2) 



kc 

where f(x) is a nondimensional function which tends to 1 as x approaches 0. Comparing 
Eqs. (1 1121) . we find that their dependence on Kolmogorov's constant is different. In Eq. (j2J), 
if we assume that E(kc) follows Kolmogorov's spectrum and K = 1.6, then the constant 
multiplying n 1/,3 fc^ 4//3 is 0.27. In contrast, in Eq. the same quantity is y/Ku* m 0.48. As 
it will be shown in the later part of the paper, the choice of constant is quite crucial in LES. 
We find that v r [kc) of Eq. (JTJ) yields better numerical results compared to v t [k\kc) of Eq. (j2J). 
We believe that the calculation of renormalized viscosity is theoretically more sound than 
the calculation of spectral eddy viscosity using EDQNM approximation, therefore, former is 
more appropriate for LES than the later. 

In this paper we perform LES of fluid turbulence using renormalized viscosity. We have 
been able to apply the same procedure to magnetohydrodynamic (MHD) turbulence also, 
except that we need two renormalized parameters: renormalized viscosity and renormalized 
resistivity. The required parameters for MHD have been recently calculated by Verma 
[13J. The LES calculations have been performed on 64 3 grid, and they have been 
compared with DNS results of 64 3 and 128 3 . As described below, the inertial range in LES 
is found to be either equal or larger than that in DNS, hence our LES model is working very 
well. 
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We solve Navier-Stokes equation in Fourier space 



Q: 



= -Mkc)k 2 u>(k) - FT (u^) - ik iP (k) (3) 

where FT stands for Fourier transform. We take v* to be equal to 0.38. 

We adopt pseudo spectral method on grid size 64 3 with dt = 10~ 4 . We apply Adam- 
Bashforth scheme to integrate the nonlinear terms, and Crank-Nicholson's scheme for the 
viscous te rm . We app.y 2/3 nue to el im inate the aliasiug e rr0r s Q. We use Fast F„u rier 
Transform developed by Frigo and Johnson 2] for our calculations. Our initial condition is 
taken to be unit energy spread out in wavenumber shells from 2 to 13 with an exponentially 
decreasing distribution. The modes in a shell have equal energy but random phases, and 
satisfy divergenceless condition. The most important ingredient in our simulation is renor- 
malized viscosity, which is computed using Eq. (1) with kc = 32. Since II changes with 
time, it is computed every 0.01 dimensionless time unit. We use dissipation rate for II. We 
carry out our simulation up to 50 time units. Our LES simulation takes approximately 60 
hours on Athlon 1.7 GHz processor. 

In Fig.^we show the energy evolution as a function of time for v* = 0.25, 0.38, 0.48. The 
E vs. t plot for all three u* are overlapping. The LES results are also compared with the 
standard pseudo-spectral DNS results performed on 64 3 (DNS64) and 128 3 (DNS128) with 
identical initial condition and v$ = 2 x 10 -4 . In DNS we apply additional hyperviscous term 
l/^ 4 u(k) with k eq = 9 to overcome aliasing errors. Clearly the energy evolution for LES 
matches quite well with DNS128, but differ significantly with DNS64. Hence, our LES on 
64 3 is able to mimic DNS of 128 3 . 

In Fig. we plot E(k)k 5 / 3 Il- 2 / 3 vs. k for DNS as well as LES. A gain the normalized 
spectrum of LES matches quite well with DNS 128 at small and intermediate wavenumbers. 
Note that 64 3 LES has much larger inertial range compared to 64 3 DNS, where it is almost 
absent. We find that the wavenumber range of inertial wavenumbers (constant with k) 
is maximum for LES with v* = 0.38; in fact wavenumber range for LES is larger than 
that for DNS128. The energy spectrum for v* = 0.25 has a hump for large wavenumbers 
(underdamped case), implying that actual v* value is higher than 0.25. The spectrum 
for v* = 0.48 shows overdamped character [16]. We have done DNS128 for some more 
parameters. The trend appears to show that v* ~ 0.38 is the most appropriate choice for 
LES. Fortunately, we obtain the above value using renormalization group calculation j^. It 



is interesting to note from Fig. ^ that the temporal evolution of energy does not clearly tell 
us which v* is the most appropriate for LES. Hence we should be careful in concluding the 
appropriateness of v* using energy evolution. The energy spectrum has more information, 
and can provide us clues on the correct choice of u*. 

From Fie. 121 we obtain the numerical value of K to be 1.7 ± 0.1; this value is close to 
the theoretically calculated va. U e !.6 QQ. Hence, the _. 1Z ed v 18 cos,ty predicted by 
Verma Q appears to be consistent and provides us a very good scheme for LES. 



and Miiller 



For LES of magnetohydrodynamic (MHD) turbulence, Agullo et al. 
and Carati [3, [l^ applied dynamic gradient-diffusion subgrid model. The forms of eddy- 
viscosity and eddy-resistivity are derived using dimensional arguments, but the constants 
are calculated using dynamical LES procedure. Their results match very well with DNS 
counterpart. In one of their main models, turbulent viscosity v t ~ /^(e^) 1 / 3 and turbulent 
resistivity r\ t ~ T 4 / 3 ^*^) 1 / 3 , where I is the resolvable length scale on the LES grid, and e K 
and e M are kinetic and magnetic energy dissipation applied by the subgrid scale respectively. 
Zhou et al. |20j have studied subgrid scale and backscatter model for MHD turbulence 

n q 

using EDQNM closure scheme. Verma pj, |2j has also calculated the above parameters 
using renormalization group procedure. Simple calculations show that turbulent dissipative 
parameters of Verma differ significantly from those of Agullo et al. and Miiller and 
Carati [rsl. figj] . as well as from those of Zhou et al. j^J. In the following discussions we will 
compare the LES results from the above three approaches. 

For MHD turbulence we apply the same LES method as described for fluid turbulence 
using renormalized parameters. The pseudo-spectral method to solve MHD equations is 
very similar to that of fluid turbulence. We also confine ourselves to zero cross helicity, i.e., 
(u • b = 0), and zero mean magnetic field. The difference of LES and DNS is in the values of 
viscosity and resistivity. In DNS we take uq = 0.00015 and 770 = 0.00015 with hyperviscosity 
and hyperesistivity parameters k eq = 7. However, in LES we take v(kc) = v r {kc)-, and 
v(kc) = T]r(kc), where kc is the cutoff wavelength. The renormalized viscosity ^(fcc), and 
renormalized resistivity rj(kc) are taken from Verma as 

u r (k c ) = (r) 1 /^ 1 / 3 ^ 3 ^ (4) 

Vr(kc) = (inv^/v/y. (5) 

Here K u is Kolmogorov's constant for MHD, n is the total energy flux, and u*, rf are renor- 
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malized parameters. The parameters ^*,?7*, and K u depend on the Alfven ratio ta, which 
is the ratio of kinetic and magnetic energy. In our decaying MHD turbulence simulation, we 
start with unit total energy and ta = 8.0. The ratio of magnetic to kinetic energy grows as 
a function of time as expected. Therefore, we need to compute the renormalized parameters 
for various values of r&. The parameters have been calculated using the procedure described 
in Verma j^j], and they are shown in Table 1. We use the appropriate v* and if given in 
the table for our simulations. The energy cascade rates are computed using Fast Fourier 
Transforms [15]. We take v r {kc) and t] r (kc) from Eqs. (@J0)J). The energy flux II changes 
with time; we compute II dynamically every 0.01 time-unit. We carried out LES for MHD 
up to 25 nondimensional time units, and it took approximately 55 hours. 

The evolution of kinetic and magnetic energies are shown in Fig. |3]as a function of time. 
The evolution of kinetic energy using LES is quite close to that using DNS. However, the 
evolution of magnetic energy does not match very well. Comparatively, LES of Agullo et 
al. l3] and Miiller and Carati l^J Q] yield a better fit to the temporal evolution of energy. 
Fig.0]shows the energy spectra of kinetic and magnetic energies for ta = 0.5 at 27 time units 
of DNS and 12 time units of LES. We find that the energy spectra calculated in LES matches 
quite well with that in DNS. The Kolmogorov's constant as indicated by the straight line in 
upper part of Fig.lUis found to be 1.8 ±0.2, which is close to the theoretical value calculated 
b 1 1 We eonCude that the LES based on _h Z ed pa—s of VerrnM is 
quite good. Our numerical results are comparable with results of Agullo et al. and 
Miiller and Carati [Lja, 1121]. However, we believe that our parameters, which are based on 
field-theoretic calculations, are on a somewhat stronger footing as compared to those used 
in earlier LES methods. 

To conclude, we have devised a LES procedure for fluid and MHD turbulence in Fourier 
space using the renormalized parameters. We take renormalized parameters from Verma Q, 
1^1 and carry out LES for 64 3 grid. When LES results are compared with DNS of size 
128 3 with the same initial conditions, we find that our LES results on energy evolution and 
spectra match quite well with the DNS results, except for the temporal evolution of magnetic 
energy. The inertial range of LES is much larger compared to DNS of the same size. Our 
results shows that substitution of renormalized parameters for eddy viscosity in LES yield 
excellent results. Hence, we demonstrate the usefulness of renormalized parameters in LES 
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calculations. 
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Figure Captions 



Fig. 1 Temporal evolution of energy in fluid turbulence using DNS and LES. The figure 
contains Energy(E) vs time plots for DNS 128 (solid line), DNS (DNS64), and three LES 
runs using v* equal to 0.38 (LES1), 0,25 (LES2) and 0.48 (LES3). The evolution in LES 
for all the three v* is quite close to DNS 128, but not to DNS64. 

Fig. 2 Energy spectrum for fluid turbulence is calculated using DNS and LES. The 
figure contains plots of normalized energy spectrum E'{k) = E(k)k 5 / 3 e~ 2 / 3 with wavenum- 
ber k after 50 time units for DNS128, DNS64, and three LES runs using v* equal to 
0.38 (LES1), 0,25 (LES2) and 0.48 (LES3). We get the best inertial range for u* = 0.38. 
The Kolmogorov's constant is found to be 1.7±0.1. DNS64 run has hardly any inertial range. 

Fig. 3 Temporal evolution of total kinetic and magnetic energy in MHD turbulence using 
DNS and LES. The kinetic energy matches quite well in both the schemes, but magnetic 
energy evolves somewhat differently. 

Fig. 4 Plots of normalized spectra E'{k) = E{k)k b l 3 e~ 2 l 3 with wavenumber k for MHD 
turbulence. The straight line shows the value of K for LES run. 
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TABLE I: The values of renormalized parameters 


for viscosity (V*) and resistivity rf in 
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FIG. 1: Temporal evolution of energy in fluid turbulence using DNS and LES. The figure contains 
Energy {E) vs time plots for DNS 128 (solid line), DNS (DNS64), and three LES runs using v* equal 
to 0.38 (LES1), 0,25 (LES2) and 0.48 (LES3). The evolution in LES for all the three u* is quite 
close to DNS128, but not to DNS64. 
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FIG. 2: Energy spectrum for fluid turbulence is calculated using DNS and LES. The figure contains 
plots of normalized energy spectrum E'(k) = E(k)k 5 ^e^ 2 ^ 3 with wavenumber k after 50 time units 
for DNS128, DNS64, and three LES runs using v* equal to 0.38 (LES1), 0,25 (LES2) and 0.48 
(LES3). We get the best inertial range for u* = 0.38. The Kolmogorov's constant is found to be 
1.7 ± 0.1. DNS64 run has hardly any inertial range. 
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FIG. 3: Temporal evolution of total kinetic and magnetic energy in MHD turbulence using DNS 
and LES. The kinetic energy matches quite well in both the schemes, but magnetic energy evolves 
somewhat differently. 
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FIG. 4: Plots of normalized spectra E'{k) = E(k)k 5 ^ 3 e 2//3 with wavenumber for MHD turbu- 
lence. The straight line shows the value of K Q for LES run. 
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